QCD corrections to J/ip production in association with a W-boson at the LHC 



Li Gang(">, Song Mao**', Zhang Ren-You**', and Ma Wen-Gan*'') 
School of Physics and Material Science, Anhui University, Hefei, Anhui 230039, P.R.China 
Department of Modern Physics, University of Science and Technology of China (USTC), 
Hefei, Anhui 230026, PR.China 

We calculate the next-to-leading order (NLO) QCD corrections to the Jltp + W production at the LHC, and 
provide the theoretical distribution of the JltJ/ transverse momentum. Our results show that the differential cross 
section at the LO is significantly enhanced by the NLO QCD corrections. We believe that the comparison 

between the theoretical predictions for the Jjij/ + W production and the experimental data at the LHC can 
provide a verification for the colour-octet mechanism of non-relativistic QCD in the description of the processes 
involving heavy quarkonium. 
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The study of heavy quarkonium is one of the most 
interesting subjects in both theoretical and experimen- 
tal physics, which offers a good ground for investigat- 
ing Quantum Chromodynamics (QCD) in both pertur- 
bative and non-perturbative regimes. The factorization 
formalism of nonrelativistic QCD (NRQCD) 1 1 ] pro- 
vides a rigorous theoretical framework to describe the 
heavy-quarkonium production and decay by separating 
the transition rate (production cross section or decay 
rate) into two parts, the short-distance part which can 
be expanded as a power series in and calculated 
perturbatively, and the long-distance matrix elements 
(LDMEs) which can be extracted from experiments. 
The importance of the LDMEs can be estimated by us- 
ing velocity scaling rules [2]. A crucial feature of the 
NRQCD is that the complete structure of the quarko- 
nium Fock space has been explicitly considered. 

By introducing the color-octet mechanism (COM), 
the NRQCD has successfully absorbed the infrared di- 
vergences in P-wave |lll|3i|4|] and D-wave [5,^] decay 
widths of heavy quarkonium, which can not be han- 
dled in the color-singlet mechanism (CSM). The COM 
can successfully reconcile the orders of magnitude of 
the discrepancies between the experimental data of 7/ ^ 
production at the Tevatron ff\ and the CSM theoretical 
predictions, even if they have been calculated up to the 
NLO. The DELPHI data also favor the NRQCD COM 
predictions for the yy — » J/ifr + X process |81,|9(]. Sim- 
ilarly the recent experimental data on the J/^ photo- 
production of HI ilOll are fairly well described by the 
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complete NLO NRQCD corrections |11], and give a 
strong support to the existence of the COM. However, 
the observed cross sections for the double charmonium 
production at B factories are much larger than the 
LO NRQCD prediction [ 1 3,1 . This discrepancy can be 
resolved by considering the CSM NLO QCD correc- 



tions 1 14] and the relativistic cou'ections |15[l without 
invoking the color-octet contributions 11611 . Further- 
more, the J /iff polarization in hadroproduction at the 
Tevatron ] 17] and photo-production at the HERA fl^ 
also conflict with the NRQCD predictions. Therefore, 
the existence of the COM is still under doubt and far 
from being proven. The further tests for the CSM and 
COM under the NRQCD in heavy quarkonium produc- 
tion are still needed. 

In order to test the COM, it is an urgent task to study 
the processes which significantly depend on the pro- 
duction mechanism. The J/iJ/ production associated 
with a W boson at the LHC, pp ^ J/iJ/ + W + X , 
can serve as a such kind of process ifioll . For this pro- 
cess, only the i color-octet (the cc[^5 j^'] Fock state) 
provides contribution at the leading-order (LO). Even 
including the NLO QCD corrections up to the alv^ or- 
der, there are only color-octets cc['5q^^], cc[-'5j^^] and 
cc[^Pf] (J = 0, 1,2), but no color -singlet contribution 
exists in the pp — > J/i// + W + X process. Therefore, 
the J/t// + W production at the LHC is an ideal ground 
to study the COM. 

As we know, the NLO QCD coiTections to quarko- 
nium production are usually significant ifl^l |20[ 21]. 
We should generally take the NLO QCD corrections 
into account in studying the COM and the universality 
of the LDMEs. In this paper, we calculate the Jjif/ + W 
production at the LHC up to the c^v^ order within the 
NRQCD framework by applying the covariant projec- 
tion method ]3], and present the theoretical prediction 
of the pj'^ distribution. The LO cross section for the 
parent process pp — > J/i// + W"^ + X involves the con- 
tributions of the following partonic processes. 



ud ^ cc[^Sf] + W^, du^ cc['S\'"] + W-. (1) 
Since the cross sections for the ud — > cc[''5j^'] + W'^ 
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FIG. 1: Representative Feynman diagrams for pp — > J lip + 



and dii — > cc[^5j^'] + partonic processes are the 
same due to the CP-conservation, we present only the 
detailed description for the calculation of the partonic 
process u(pi)d(p2) — > cc\^S^f^]{pi) + W^(p^). The 
tree-level diagrams for this partonic process are drawn 
inFigs[na)-(b). 

In the nonrelativistic limit, by applying the covariant 
projection method we obtain the differential cross 
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section for ud — > cc[ 5 j ] + expressed as 



d&o 



{f + M^) + iu + f + 
+ 2s{t + u)\}, (2) 
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where s = {pi+p2) , f = (pi - P3) and u = (pi -774) ■ 
NcoiCsf) and Npoii^Sf^) refer to the color and po- 
larization degrees of freedom of cc[^5j^'] |4i], respec- 
tively. Then the LO cross section for the pp J/t// + 
-I- X process is 



r(0) 



= J dxidx2d&o[Gu/A(xuHf)G^/B{x2,pf) 
+ (A^B)], (s^xiX2s), (3) 



where G„ j/^ g are the parton distribution functions 
(PDFs), and A, B represent the two incoming protons 
at the LHC. 

In calculating the NLO QCD corrections to the 
pp — » J/t// + + X process, we should consider 
both the virtual correction and the real gluon/light- 
quark emission correction. The virtual corrections 
only come from the cc[^5j^^] Fock state, while the 
real gluon/light-quark emission correction involves the 



contributions of the cc[^5q^''], cc[-^5j**'] and cc[^pj^] 
(y = 0, 1, 2) Fock states. In our calculations, we adopt 
the dimensional regularization (DR) scheme to regular- 
ize the UV and IR divergences, and the modified min- 
imal subtraction (MS) and on-mass-shell schemes to 
renormalize the strong coupling constant and the quark 
wave functions, respectively. 

There are 41 virtual QCD one-loop diagrams for 
the subprocess — > J/i/> + W'^, which include self- 
energy (12), vertex (10), box (7), pentagon (2) and 
counterterm (10) diagrams. We present part of these 
diagrams in FigslTJc)-(l). There exist ultraviolet (UV), 
Coulomb and soft/collinear infrared (IR) singularities 
in the virtual correction. The UV singularities are can- 
celed by the counterterms of the strong coupling con- 
stant and the quark wave functions after the renormal- 
ization procedure. But the QCD one-loop amplitude 
of the partonic process ud cc[^Sf ^] + still con- 
tains Coulomb and soft/collinear IR singularities. The 
IR and Coulomb singularities in the virtual correction 
can be expressed as 
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where 

Al = -2Cf, 
Cf = 4/3. 
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The soft/collinear IR singularities can be canceled 
by adding the contributions of the real gluon and 
light-quark emission partonic processes, and redefin- 
ing the parton distribution functions at the NLO. For 
the Coulomb singularities, they can be canceled af- 
ter taking into account the corresponding corrections 
to the operator < 0^^'^\^Sf\ >. We use the expres- 
sions in Refs. [22- 24] to implement the numerical eval- 
uations of IR-safe one-point, 2-point, 3-point, 4-point 
and 5-point integrals. In our calculations, only two dia- 
grams, FigsUlf) and (h), contain the Coulomb singular- 
ities, which are regularized by a small relative velocity 
V between c and c [25]. We adopt the expressions in 
Ref.[26] to deal with the IR-divergent Feynman inte- 
gral functions. 

The real gluon emission process provides three types 
of corrections, which correspond to the contributions 
from the cc[i5[,^^], cc[^Sf] and cc[^pf] (J = 0, 1,2) 
Fock states, respectively. The real gluon emission cor- 
rection to the ud — > cc['5g^']-i-W^ subprocess is free of 
divergence, and can be numerically calculated by using 
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the Monte Carlo method. For the cc[^S f ^] + +g pro- 
duction, it contains both soft and colUnear IR singulari- 
ties which can be conveniently isolated by adopting the 
two cutoff phase space slicing (TCPSS) method [27]. 

In adopting the TCPSS method, we should introduce 
two arbitrary cutoffs, 6s and 6c- The phase space of 
the ud cc[^Sf ^] + + g subprocess can be split 
up into two regions, £5 < 6s (soft gluon region) 
and Es > 5.5 VI/2 (hard gluon region) by soft cutoff 
^5. Furthermore, the hard gluon region is separated as 
hard collinear (HC) and hard non-collinear (HC) re- 
gions. The HC region is the phase space where -fis 
(or -t2i)< 6cS (fi5 = ipi - pif and = (pi - psf)- 
Therefore, the cross section for the real gluon emission 
subprocess can be expressed as 

^(^5f) = 



^(^5f). 
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The cross section for the subprocess ud — » cc[-^5 -H 
+ g in the hard non-collinear (HC) region is free of 
divergence, and can be numerically calculated by using 
the Monte Carlo method. The differential cross section 
for the subprocess ud ■ 
region can be expressed as 



cc[^5f' ] + +gm the soft 
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where g(pi,Pj) are soft integral functions defined as 



For the subprocess ud — > cc[-'5j^'] + W'^ + g only 
the gluon radiation from initial particles can induce 
collinear singularities. The differential cross section, 
d&^'^, can be written as 



do-: 
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where Pu„(z, e) and P^^{z, e) are the D-dimensional un- 
regulated (z < 1) splitting functions related to the usual 
AltarelU-Parisi sphtting kernels Isill . Pii(z,e)(i - u,d) 
can be written explicitly as 



Pn(z,e) = 

Pii(z) = 
P'n(z) = 



Puiz) + eP'iiiz) 
-Cf{1-z) 



(i — u, d). 
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cc[^P'^'] + W^+g, 



As for the partonic process ud — > ccl^Pf*] 
it contains only the soft singularities. Using the TCPSS 
method mentioned above, we split the phase space up 
into soft gluon region and hard gluon region by adopt- 
ing the cutoff 6s. Then the cross section for the partonic 
process ud — > cc\^P^^^] + + g can be expressed as 
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The cross section &^(^pj^) is finite and can be evalu- 
ated in four dimensions by using Monte Carlo method. 
The differential cross section in the soft region for the 



P5 



^(Pi ■ P j) process ud 



cc[^Pf] + 



W+ +g, d&l(^Pf), can be 



_(Pi ■ Ps){Pj ■ P5) expressed as 
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Then we can get 
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• 4-mt/Ei and £3 = — . 

The real light-quark corrections to the subprocess 
ud ^ J/i// + arise from the partonic processes 

g{pi)u{d){p2) cc[n](pi) + W^ip4) + d{u){p5),{l5) 
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where n = f \ and "'Py Some of the Feynman 
diagrams for these partonic processes are presented in 

Figsino)-(p). 

The real Hght-quark partonic processes with 
n = ^Sf-' and ^Pf' contain no singularities, so we 
can perform their phase space integrations by using the 
general Monte Carlo method. The partonic processes 
gu{d) — > cc[-'5j^*] + + d(u) contain only the ini- 
tial state collinear singularities. By using the TCPSS 
method, we split the phase space up into two regions, 
collinear region and non-collinear region. 



Af (a ^ ag) 
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+ A^"(a ^ ag) 



where 



^Gj/g(x2,fif) + (A «-» B)\dx\dx2, 



Af (a ^ ag) = Cf(21n5, + 3/2), 
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djCsf) = a-':esf^) + ^U'sf^) (q = u,d). (16) and 



The cross section in non-collinear region, &^Csf) IS 
finite and can be evaluated in four dimensions by using 
the Monte Carlo method. The differential cross sec 
tions for the subprocesses gu(d) 
d(u) can be written as 



s 



Ga/H(x,IJ.f) = 2j I —Ga'/H(x/y,tJ.f)Paa'(y), 
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X Sc'[Pd(u)g(z, e)Gg/A(xi/z, /i/)G„(rf)/B(x2, fif) 
dx\dx2 (17) 
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We can see that the summation of the soft (Eq.®), 
collinear (Eq.(l20t). and UV renormalized virtual cor- 
rections (Eq.®) to the pp J/t// + + X process, 
dcrl + dcr^"" + dcr^, is soft/collinear IR- and UV-finite, 



with q - u,d and Pqg(z, e) can be expressed explicitly i-e., 

as 



0, 



P,g{z,e)^P,g{z) + eP'^^{z). 

P,g{z)^^-[z^ + {\-z)\ 

P;/z) = -z(l-z). (18) 

To obtain an IR-safe cross section for the pp — » 
y/i^ + +XvipX.o the NLO, we should take into ac- 
count both the NLO QCD counterterms of the PDFs 
and the NLO QCD corrections to the operator < 
^y/^[3^(8)j ^ rpj^g counterterms of the PDFs 



^^2 

A\ + Af (m ^ ug) + Af (J ^ dg) = 0. 



(24) 



+ do-^, still 
the 



are expressed as [3. 



6Gi/A(x,pf) = - 



as r(l - e) 
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X j jPu(z)Gj/A(x/z,fXf). (19) 

By adding the contributions of the PDF counterterms 
and the real gluon/light-quark emission collinear cor- 
rections shown in Eqs.lfTTI) and ( fTTb .we obtain 



However the above result, dcr^^ + dcr^"" 
contains the Coulomb singularity. Furthermore 
corrections contributed by the ccl^Pf*] (J - 0, 1,2) 
Fock states, d&^i^Pf^), contain soft IR singularities 
too. We can eliminate these singularities by taking 
into account the NLO QCD corrections to the opera- 
tor < O-'^'^'l^Sf^] >. In this paper we use the method 
in Ref.||4ji] to deal with these singularities. In Fig|2]we 
present the IR and Coulomb singularities structure in 
the NLO QCD calculations for the pp J/i// + W^ +X 
process. We have checked analytically that all the IR 
and Coulomb singularities are canceled in the final re- 
sult. 

The final result for the process pp — > + W^X up 
to the NLO consists of three parts of contributions. 



O" rotal - (Ti^m H- (T, ^(8) + O", 
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where cr3„(g) can be divided into two parts: a two-body 
term cr'^'g, and a three-body term o-f\g,. The two-body 
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term, cr is expressed as 



0-^2) = o-(0) + 



Us f 

— I dxidx2dd-o{Gu/A(xufif) 

X G3/B{x2,iJf)[A^o+Al+Al^(u^ug) 
+ A^id^ dg)] + GuiA{xutif)GjiB{x2,iif) 
+ G„/A(jci,yU/)Q/_Bfe,yU/) + (A «^ B)). (26) 



And the three-body term, cr^^\^^, is written as 



,(3) 



^HC^3c(8)^ 



(27) 



0"3„(8) also can be divided into two parts: a two-body 



term cr*^' s) and a three-body term cr\^\i^ ■ The two-body 



,(3) 
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term, cr'^'g,, is expressed as 



3 p(8> 



J" dxidx2d&g (^Pf^){Gu/Aixi , n f)Gj/g(x2, /J. /) 



H- (A B)]. 
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And the three-body term, cr\ is written as 



(28) 



(29) 



Finally, after taking into account the NRQCD NLO 
corrections to the operator < 0''^'^[^Sf ^] >, the con- 
tributions of the and ^Pf' states are finite. As for 
the contribution from ' S f ^ state, it contains no singu- 
larity and only involves a three-body term crf \g,, which 
can be expressed as 



(30) 



As a check of the correctness of our calculations, the 
independence of the cross section part of o-,^m -i-cr,^(8) -h 
cr,^(8i, on the two arbitrary cutoffs, i5.5 and 6c, has been 
numerically verified. 

For the pp ^ J/t// + W + X process at the LHC, 
we take CTEQ6L1 PDFs with an one-loop running 
as in the LO calculations, and CTEQ6M PDFs with 
a two-loop cKj in the NLO calculations |33]. For the 
QCD parameters we take the number of active flavor 
a&tif - 3, and input Aq^^ - 249 MeV for the LO and 



A' 



(3) 
QCD 



389 MeV for the NLO calculations lUJl, re- 



spectively. The masses of the external particles and the 
fine structure constant are taken as niw = 80.398 GeV, 
mc = mj/^/2 = 1.5 GeV and a = 1/137.036. The 
sine squared of the Weinberg angle is expressed as 
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FIG. 2: The IR and Coulomb singularities structure in the 
NLO QCD calculations for the pp — » J/if/ +W'^+X process. 



= 1 - m^/m|, where mz - 91.1876 GeV. The 
renormalization, factorization and NRQCD scales are 
chosen as yUr = A*/ = mj and //a - iric, respectively. 



where mr - 



+ nt^ji^ is the Jjij/ transverse 



mass. Since the < 0^^'l'[^Pf] > (7 = 0, 1,2) LDMEs 
satisfy the multiplicity relations 

< O-'I'^i^Pf] (27 -H 1) < 0^^'^i^P^^^] >, 



we adopt the LDME < O^I'I'i^Sf] 
10"^ GeV^ and the linear combination 



Mi 



> +- 



>= 2.73 X 



wifli m!'''' = 5.72 X 10-' GeV" and r = 3.54 as 
the input parameters, which were fitted to the Teva- 
tron RUN-I data by using the CTEQ4 PDFs and have 
taken into account the dominant higher-order eff'ects 
due to the multiple-gluon radiation in the inclusive J/ tfr 
hadroproduction |34]. Then the < 0-"'''[^Sf^] > and 
< 0-"'''pPj^^] > LDMEs are fixed by the democratic 
choice I 35I1 

mi 2 

In the calculations of the real corrections, the two 
phase space cutoff's are chosen as Ss = 10"^ and 6c = 
6s/50, and the invariance of the 6s value running in the 
range of 10"^ - 10"^ is checked within the error toler- 
ance. Considering the validity of the NRQCD and per- 
turbation method, we restrict our results in the range of 
p'^"^ > 3 GeV and \yj/^\ < 3. 
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FIG. 3: The LO and NLO QCD corrected distributions of 



p'j'^ for the pp J 14'+ M'* +X process, and the contributions 
of the cc ^Sf^ and cc ^Sf\ Fock states up to NLO at the 



LHC. 



In Fig 01 we present the LO and NLO QCD cor- 
rected distributions of pj'^ for the pp Jlij/ + +X 
process at the LHC. For comparison, we also depict 
the contributions of the cc['5q^'] and cc[^5[*'] Fock 
states in the figure, while the total contribution of the 
cc[-'P®] {J - 0, 1,2) Fock states is negative and will 
be drawn in Fig|4] From the Figl3] we can see that 
the differential cross section at the LO is significantly 
enhanced by the QCD corrections. In the range of 
3 GeV < p-'J'^ < 50 GeV, the /T-factor, defined as 



K = 



da-" 



-rwl^, is in the range of [3.09, 4.31], and 

reaches its maximum when p^J.''' = 3 GeV. In Figg] we 
present the various contributions of the cc[^Py '] Fock 
states. The solid, dashed and dash-dotted curves repre- 
sent the total contribution of the cc\^P^j^] Fock states, 
the real gluon emission correction to the cc\^P^j^] + W 
production and the real light-quark emission correction 
to the cc\^P'"f] + W production, respectively. In Figs|3] 
andlH we can see that there exists a compensation be- 
tween the contributions of the S and P states. The con- 
tributions of the cc['5[)**'], cc[''5 j^'] Fock states and the 
real light-quark emission correction to the cc\^P^j^] + W 
production are always positive, while the real gluon 
emission correction to the cc[^Py '] + W production is 
negative. 

In above calculations, we mainly consider the di- 
rect production of y/i/r mesons up to the NLO. But 



FIG. 4: The ratio R versus Pj 
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the y/i/r also can be produced indirectly via radiative or 
hadronic decays of heavier charmonia, such as Xcj and 
1^' mesons. The respective decay branching fractions 
areBOrcO ^ ^M+r) = (1.28±0.11)%,B(;r.i ^ + 
r) = (36.0+1.9)%, Bixci J/i/^ + y) = (20.0+ 1.0)% 
and B(^' J/if/ + X) = (57.4 ± 0.9)% [36]. The cross 
sections for these four indirect production channels can 
be obtained approximately by multiplying the direct- 
production cross sections for the respective intermedi- 
ate charmonia by their decay branching fractions. At 
the LO, only the cc[''5j^*] color-octet contributes to 
the productions of the charmonia associated with a W 
boson at the LHC. With the multiplicity relations of 
LDMEs < CF^'[^sf '] >, < 0^"[^Sf ^] >= (27 + 1) < 
qxm cross sections for the productions of 

the charmonia associated with W can be expressed as 



d&(pp H + W) - d&Q 



(31) 



where H represents the intermediate charmonia i/r' or 
Xcj- In numerical calculations, we adopt the LDMEs 
< a<^<>[^f] (6.81 + 1.75) X lO-'* GeV^LMJ and 
:(8)t 



< O^'i^Sf] (2.0 + 0.2) X 10-3 GeV^CT. We find 
that the cross section for the indirect production 
is almost the same as the LO cross section for the di- 
rect production, cr'"<''"''='=' = 0.94 x cr(°>. The com- 
plete indirect y/i/r production in association with a W 
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gauge boson at the NLO will be calculated in our fur- 
ther work. 

In this paper we investigate the NLO QCD correc- 
tions to the Jl^+W production at the LHC, and present 
the numerical results of the differential cross section 
of the Pj'^ and the contributions of the different Fock 
states for the pp J/iff + W+X process up to the QCD 
NLO. We find that the differential cross section at the 
LO is significantly enhanced by the QCD corrections. 
The numerical results show that there exists a negative 
real gluon emission correction to the cc[^Pf^] + W pro- 
duction. The LO differential cross section for the pro- 
cess pp ^ J /if/ + W- + Z in the low p^^'^ region is 
heavily enhanced, and the NLO QCD corrected differ- 
ential cross section can reach 0.49 pb in the vicinity of 
p^J'^ ~ 3 GeV. Although the J /if/ events are difficult to 
be accepted at low pr region, we can find that the cross 
section for J/i// + W direct production at the NLO is 
about 0.81 ph with * > lOGeV . To obtain the cross 
section for the process pp — > J/i// + W + X with the 
pure leptonic decays, we multiply the cross section for 
the direct production by the branching fractions 12% 
for J/i// -» and 22% for W -» Iv. Given the in- 
tegrated luminosity of 300 fb~^ at the LHC, we could 
obtain about 6400 events. If we include the indirect 
contribution of the radiative or hadronic decays of XcJ 
and i//' mesons to J/i//, more events could be detected. 
Even taking into account the detector acceptance and 
efficiency, there are still enough J/i/z+W events which 
can be detected. We conclude that the LHC has the po- 
tential to detect the J/ifr + W production. If the J/ifi+W 
production is really detected, it would be a solid basis 
for testing the colour-octet mechanism of the NRQCD. 
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